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ABSTRACT 


‘> A method is derived for finding the linear response and 
loading transfer functions for the lateral aerodynamic case 
of airship flight through atmospheric turbulence. The func~ 
tions obtained are in a form that can be applied to the 
various spectral analysis methods used to predict surviva~ 
bility currently employed by designers. A numerical example 
uging the USS AKRON (ZR=4) is presented. The results show 
that peak motion response and loading occur when the an- 
countered spectral component has a wavelength equal to the 


airship length, and that simple feedback of heading angle 


does not significantly decrease this peak. 
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BM 


NOMENCLATURE 


hull cross-sectional area 


total buoyancy force of the airship, 
p*g:voluma 


hull banding moment about the vartical axis 
np eae along the hull longitudinal 
axis 


longitudinal characteristic length of the 
airship 


mean chord of fin 


matrix of coefficients from equation 51 
nondimensional aerodynamic force in the y~- 


direction, Cy = 2¥/oU.,?s 


nondimensional aerodynamic rolling and yawing 
moments respectively, (C,, Cy) - 2(L,N) /pU,286 


nondimangional aerodynamic Lift 
Cy, for the fins alone, no hull interference 
nondimensional shear force c. 28 (1) /oU,"8 


nondimensional banding and twisting moment 

) sawn 
respectively, (Cows Coy) ™ 2(BM,TM)/pU,"Se 
nondimensional time derivative, (¢/2U,)d( )/dt 
gravitational acceleration 


turbulence forcing functions (aquation 59) 


body~fixed coordinate measured normal to the 
hull centerline (positive up) 


: h h location of the vehicle's mass center 


om 
3 (Nom) s h location of the empennage~-assembly's mass 
center 
H(k) Sears' function corrected for finite aspect 
ratios by Filotas [Ref. 16] 
Log? Tee moments of inertia about the x- and z-axes 
I. product of inertia w.r.t. x= and z~axes 
ds doe 1. nondimensional moments and product of inertia 
=~3 
dee = 8I,./pse 
4 imaginary operator, i= /-T 
K hull potential cross-flow factor from Jones 
and DeLaurier [Ref. 6] 
Ry Ko axial and traverse apparent-mass coefficients 
Ky control gain of fin normal force to vehicle 
. azimuth angla 
1; axlally~aligned body=fixed coordinate origi- 
nating at the nose 
{ 
| 1h axial location of the buoyancy center 
| Ly axial location of the hull-fin intersection 
point 
| le axial location of the fin's aerodynamic center 
: Lia axial mass-center location of the entire 
i vehicle , 
4 L turbulent scale length 
| L rolling moment - 
| m mass of the entire vehicle, including internal 
air and gas 
ms mass of the empennage assembly 
| N yawing moment 
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HD ww 


ty> 


GREEK SYMBOLS 
to 
8 


vehicle angular velocity about the x-axis 


nondimansional value of p, 6 = (€/2U,)p 


nondimensional maximum value of 6 

vehicle angular velccity about the z-axis 
nondimensional value of r, 2 = (c/2U,)r 
nondimensional maximum value of ? 
auto-correlation function 
reference area of airship, (volume) 2/3 


hull shear force, normal to the centerline 
stabilizer reference area (planform area) 
time 

twisting moment, about the hull axis 
reference flight speed 


perturbation velocity of the vehicle's mass 
center in the y-direction 


nondimensional value of v, ? = v/U, 
nondimensional maximum magnitude of ¢ 


horizontal velocity of the atmospheric turbu- 
lence 


body~fixed wind-aligned stability axes (x 
positive forward, y positive right, z positive 
down) 

axes fixed in inertial space 


force in the y-direction 


reference aerodynamic relative angle of attack 


aerodynamic sideslip angle 


Ll 


Y horizontal nondimensional velocity of the 


spectral component - 
P maximum value of 7 | 
Ns stabilizer efficiency factor, from Jones and 
DeLaurier [Ref. 6] 
| u nondimensional mass, yu = 2m/oSe 
& axial coordinate, measured from the nose 
fe) atmospheric density 
a turbulence intensity 
es nondimensional stability root 
> roll angle 
9 maximum value of 6 
944 (2) power-spectral function for turbulence com~ 


onent v 
P g 


v azimuth angle 4 
y maximum value of a 
Rf 
| A) spectral component frequency 7 
| ry) wave number 
«| 
| SUBSCRIPTS 
| ( leaee aerodynamic force or moment terms 
‘4 (ep buoyancy terms 
Se control terms 
| : ( en mass center terms 
( emp term for entire empennage assembly 
| ( Vg atmospheric-turbulence term 
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or hull term }« 


(in inertial term 
( 3 reference equilibrium value 
‘ (5 derivative w.r.t. p 
: ( Vs derivative w.r.t. Dp 
i 
: ee derivative w.r.t. 2 
‘c (ove derivative w.r.t. Dr 
i. ie fin term 
(Vn thruster-rotor term . 
Che derivative w.r.t. v 
ie ee derivative w.r.t. Dv 
( )g derivative w.r.t. 8 
(a derivative w.r.t. Df 
SUPERSCRIPTS 
*) nondimensional term 
) derivative w.r.t. time 
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I. INTRODUCTION 


When a close look is taken of the history of airship 
flight, it becomes evident that turbulence is a prime cause 
for concern in the design of lighter-than-air (LTA) craft. 
The spectacular crashes of the airships Shennandoah, Akron, 
and Macon are perhaps the most obvious remindars of the power- 
ful effect the wind can have on these fragile vehicles. It 
is imperative that proper consideration be given to the gust 
response of an airship, both dynamically and structurally. 

At the time the great rigids were built, designers had 
only a cursory knowledge of turbulence and how to deal with 
it. Burgess (Ref. 1], for example, dedicates only one page 
to the subject, and that is mainly a warning that gusts en- 
countered in flight can place larger loads on an airship 
than any maneuver of which it is capable. He suggests using 
the standard gust analysis technique of the day--a "fixed- 
in-space" flight through a ramp shaped gust. This was the 
method used to evaluate all aircraft at the time, and pro- 
vided some measure of confidence that the structure would 
withstand the stresses of flight. Its principal weakness 
was that it took into account only the direct gust loads. 

Since that time, the methods used to analyze response to 
turbulence have been greatly improved. Over time, the ramp 


shaped gust gave way to other gust shapes, finally settling 
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on the (l-cosine) gust as standard (Ref. 2]. By the mid- 
fifties anough data had been gathered on the nature of the 
turbulent wind to derive some basic statistical relations. 
This, in turn, allowed the development of the power spectral 


method for aircraft load analysis (see, for example, Press 


i and Meadows [Ref. 3]). The spectral method has been applied 
Ht to many types of aircraft over the last twenty years. At 

: | present, both methods are used in analysis of flight struc- 
a tures, the one giving the more conservative structure de- 
termining the final design. 

While the development of turbulence modeling techniques 
ie progressed, applications to airship technology were slower 
in coming. Because the success of heavier-than-air (HTA) 
craft made the slower airship aconomically less attractive, 
LTA research lagged. In the period between World War II and 
the Arab Oil Embargo of 1973, the only significant contribu~ 
tion to the study of airships i: ‘urbulence was by Calligeros 
and McDavitt [Ref. 4]. This paper presented a method of 


analysis that allowed a stable airship to respond dynamically 


to both sinusoidal and (l-cosine) gusts. Thus, the inertial 


a ae err 


| 
| 
{ 
A and aerodynamic reaction forces of gust encounter could be 
{ 


Mal e 
ie included in the model. Further, by using the sinusoidal 

* | representation for gust shape, it was possible to apply spec- 
is Ss 

uy 


tral methods to the analysis. 
i The recent interest in LTA brought about by rising fuel 


‘a | costs has increased the research in all aspects of airship 
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flight. The spectre of the great rigids disintegrating in 
turbulence makes it imperative that the designer have ade- 
quate means to predict an airship's response to gust pene~ 
tration. Current research is aimed at supplying that means. 
DeLaurier and Hui [Ref. 5] refined the technique of Calli- 
geros and McDavitt [Ref. 4], by introducing refinements to 
the aerodynamic cross-flow model (Ref. 6] and allowing for 
stability augmentation through pilot control input. The 
model allows statistical prediction of an airship's dynamic 
response and operational lifetime for various combinations 
of speed, altitude and control gain. 

DeLaurier and Hui's work (as well as most others dealing 
with this subject) concentrated only on the longitudinal 
aerodynamic case. This thesis proposes to apply their model 
to the lateral case, enabling the response to side--force to 
be calculated. Bending and shear in the horizontal plane, 


as well as twisting moment, can then be taken into account 


when predicting airship life expectancies. 


Te RF 


ein LSS 


ae ee 


II. THE TURBULENT WIND 


The motion of the atmosphere is vary complex. Shearing 
stress between layers of differant speeds and at the ground, 
thermals caused by solar heating, weather fronts, vortex 
shedding behind obstructions and aircraft, plus many other 
phenomena, all contribute to a velocity field that is most 
d‘.fficult to describe. Of the methods available, that chosen 
will depend on the purpose for which it is used. A goal in 
the dasgign of any flight vehicle is safety, with performance 
adequate for the mission. This dictates using models giving 
reasonable estimates for the design parameters, and it may 
or may not be necessary to closely match the physical reality 
to do this. In any casa, the designer must be familiar with 
the advantages and limitations of methods available in order 
to choose wisely. What follows is a brief review of some of 
the turbulence models currently in use and how they apply to 


airship analysis. 


A. THE DISCRETE GUST 

As mentioned in the introduction, the discrete gust model 
has bean used to analyze aircraft for many years. It is 
especially good when response to the passage through a steady 


velocity gradient, such as a thermal, mountain updraft, or 


jet stream, is desired. The method has been improved 


steadily until, as Etkin [Ref. 7] points out, it has attained 
a high degree of sophistication. 

Figure 1 shows the shape for a (l-cosine) gust, where Wn 
and qn are maximum gust velocity and distance along the 
flight path of this maximum. By varying these parameters, 
the gust severity can be controlled. The value for dn is 


prescribed by the Federal Aviation Administration (FAA) as 


The size of Wa is dependent on airspeed and altitude, and is 
shown in figure 2 for three values of equivalent airspeed. 
The 25¢ wavelength was chosen because it historically couples 
with the short period pitch mode of a rigid aircraft to pro=- 
duce the largest load factors. Calligeros and McDavitt (Ref. 
4] showed that, for airships, the maximum loads occur when 
the wavelength is equal to the airship length. 

The British dictate (ARB CAR CH D3-3) that the gust para- 
meters be chosen to produce the peak response with aircraft 
flexibility taken into account. In this way, the model is 


"tuned" to the aircraft, thus assuring a conservative design. 


B. RANDOM TURBULENCE 

Extensive measurements of the atmospheric velocity field 
have been made, and the techniques involved are well estab-~- 
lished and reliable [Raf. 2]. They show that the velocity 


vector is best characterized by a random function of space 
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and time that is, in general, non~homogeneous, non- 
stationary, and anisotropic. The axact function has not 
been daveloped due to its tremendous complexity. Until 
enough data is collectad (if evar) to allow the precise 
formulation, cartain simplifying assumptions must be made 
to enable flight vehicle analysis. 

One assumption that applies everywhere except in the 
planetary bouxuié7y layer (below about 1000 ft), is that the 
| turbulence is ‘homogeneous’, that is, the statistics of the 
field do not vary through space. In the boundary layer, 
scale length and intanaity are homogeneous in the horizontal 
plane, but not vertically. Another assumption made is that 


the turbulence is 'stationary', or statistically time con- 


stant. Over the time periods of interest to flight this 
approximation is quite adequate. Also, the turbulence is 


assumed to be isotropic (again, axcept in the boundary 


layer), making the statistics invariant with orientation. 


eo 


One last simplification used to model atmospheric tur- 


bulence is the ‘frozen field' or Taylor's hypothesis [Ref. 


8]. The change in the velocity field perceived by an air- 
a craft as it passes with speed Us through the air is given by 


i the substantial derivative 


ce 3 Bed = Hed + ue He? | 


SS 


A Pe 


Taylor's hypothesis states that for all but the smallest 
i values of U,, the sacond term dominates and the first may 
be ignored. The result is that the correlations and spectra 
reduce to threa-dimensional functions of space only. Phy~ 
sically, this means the velocity field is ‘frozen' in time, 
and the changes are due only to displacement. 
In an effort to specify an acceptable lower limit on Us? 
for which the frozen field applies, Dobrolenskiy [Ref. 9] 
cites studies comparing records of turbulence spectra gath~ 
erad by captive balloon and an aircraft flying nearby at the 
same time. Within the margin of error, the two are statis-~ 
tically quite similar, and he concludes the lower limit on 
U5 is comparable to the convection velocity (for all prac- 
tical purposes, the mean wind speed). Etkin [Ref. 7] points 
out that the vehicle speed can be as low as one-third the 
wind speed for good results. Note that the only vehicles 
| capable of less than this velocity are LTA and VTOL craft, 
and then only when they are convected downwind with the air- 
mass. In hover, or upwind flight, the hypothesis holds. 
For that small portion of the flight envelope in which it 


{ 
| 
| 
t 
aoa 
4 do@s not, the forces generated on an airship's structure are 
small and, therefore, do not present a problem. 


The techniques for dealing with isotropic, frowan turbu- 


lence are well known [Ref. 8]. For those not familiar with 


the mathematical background necessary to deal with the 


ZL 


subject in depth, Chapters 2, 3,and 13 of Etkin [Ref. 10] 
will provide a good primer. 
1, Turbulence Organization 

Anyone who has seen leaves swirl on an autumn day, 
or watched someone blow smoke rings, has an intuitive under~ 
standing that the motion of the atmosphere is not completely 
arbitrary. What happens at one location effects conditions 
at another. Fuidynamicists characterize this interdependence 
by using expressions relating the stress and strain in the 
fluid. The technique has been applied to turbulence (Ref. 
11] with some success in predicting the actual velocity field 
of a boundary layer type flow. 

For flight vehicle analysis, a more convenient method q 
of specifying the velocity field is the spactral decomposi- 
tion of the three-dimensional homogeneous vector field (Ref. 
12). Figure 3 shows an aircraft flying through a (two- 
dimensional) sinusoidal wave of shearing motion. The velo- 
city change from the mean is given (for the lateral gust 


component) by 


L(Q,x' + Noy") 
dv (x",y") “ @ de, (1) 


where Q, and Qo are the wave number components in the x' and 
y' directions respectively, and Co ds the complex amplitude 


of the lateral component, If the vehicle penetrates the 
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field with velocity Use the coordinatas become x' = x + Ust. 
and y' = y in the body fixed system. Equation (1) then 


becomes 


iN,WU.t L(QSK + Qyy) 
dv (x,y) a @ vo e 1 2 de. (2) 


The air velocity over the vehicie is then periodic with wave- 
lengths (20/24 9) and frequencies (25 2U,/2). The total 
field is mada up of the superposition of these spectral 
components, much as a Fourier series represents a random 


scalar. 


2. Probability Distribution and Spectra 


With the expression for a single spectral component 
available, the next difficulty is determining the probability 
distribution of the individual frequencies. The power spen- 
tral density of a time varying function, X(t), is defined 


(in terms of wave number) as 


1 


y 
To i X(t,2, ARs dt (3) 


p(2) = Lim 
AN~+0 


Ton 


where $(2) is expressed in (ft/sec) */ (radians/#£t) , and T is 
the duration over which X(t) is measured. The value is usu- 
ally computed by taking the autocorrelation function R(t) 
(Ref. 10: Chapter 2], and performing a Fourier transforma~ 


tion, thus 
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1, 


v(ay = ge | Rinve! ar (4) 


~00 
Figure 4 [Ref. 2] shows the spectra of three different 
meteorological conditions. While varying in detail, they 
all exhibit the same decreasing trend at higher frequency. 
The vertical dimension is a measure of the intensity of 

the turbulence at that particular frequency, and the square 
root of the area under the curve a measure of the overall 
rms gust velocity [Ref. 13]. As a point of practicality, 
only the area under the actual measured curve is included, 
because values of the spectrum at higher frequencies con- 
tribute Little to the response of an aircraft, 

From analyzing larga numbers of samples, it is ap- 
parent the probability distribution is non-Gaussian [Ref. 
14], with very high and very low values of intensity occur- 
ring more frequently than predicted by a normal distribution. 
However, the vast majority of values do fall on a Gaussian 
curve, sO it is reasonable to use this assumption for most 
applications. This is very beneficial because, whereas 
Gaussian input to a linear system results in Gaussian out- 
put, response to non-Gaussian input is, in general, unknown. 
In the analysis of flight vehicles, linear system models are 
used axtensively, making the assumption of normal distribu-~ 
tion most desirable. In order to account for the large 


gusts omitted by this type model, the (l-cosine) method can 
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be employed. This is the practice recommended by many cer- | 
tifying agencies. 


Two Gaussian models in current use are the Dryden 


spectrum 
I ( 
2s ~2.2 
L 1+ 3L°9 
boo (Q) = St (5) 
a7 SG eae y2 


and the von Karman spectrum 


2-1 + 8 (1.339%) 2 
oa L 3 
Sag (2) SS - (6) 
[1 + (1.339E9) 2} 20/20 


The first was dyevaloped to define turbulence spectra in wind 


tunnels. It is the simpler of the two, but not as accurate 


as in the past, but is mentioned here due to its historical 
\ significance and the large number of references to it in the 


literature. 


| 
| 
| 
| 
| 
as the second. For that reason it is not used as much today 
| 
Today, the von Karman spectrum is used almost uni- | 
versally. In aquation (6), o is the rms turbulence intensity, | 
E 

ss 2 

and L is the “scale length"--a measure of the average eddy | 


size encountered. Testing has shown that the model is a : 


\ reasonable fit to all levels of turbulence. Figure 5 is a | 
i 


cr Be Ra err ST 


plot of one set of experimental data along with the predicted 
turbulence spectra for severe storm conditions using various 


' 

| 

: 

values for L. As can be seen, the agreement is quite good. | 
| 
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The values for o and L are variable and appear to be 
functions ot altitude. In addition, two standard categories 
of intensity are defined--"storm" and "non-storm". Table I 
[Ref. 15] lists values for non-storm (b,) and storm (b,) 
intensities, as well as scale length, currently used by NASA 
for horizortal atmospheric flight. In this table, the values 
Pp, and po are the probabilities of encountering non-storm 
and storm turbulence, respectively, at the altitude speci- 
fied, Note that the values of L given for flight below 1000 

* are representative values that are probably low, giving 
conservative (high) numbers of load exceedancas per unit 
length of flight. These values are usable for structural 
analysis, but inappropriate for control system studies of 
flight simulation where the vertical inhomogeneity must be 
taken into account. 

Equation (6) is the expression for the transverse 
Spectrum needed to analyze the longitudinal response to 
vertical gusts. When dealing with lateral aerodynamics, such 
as the analysis done in Chapter III, it is necessary to use 


the longitudinal spectrum 


22 2 
w [lL + (1.339L2) *] 


where o and i are the game as in the transverse case. The 
term "longitudinal" refers to the variation in gust velocity 


parallel to the direction of the mean wind, whereas 
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"transverse" refers to the perpendicular direction (both 


vertically and horizontally), thus 


It should be obvious from this discussion that flight direc- 


tion relative to the mean wind becomes important when dealing 


with lateral aerodynamics. 
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IIIT. METHOD OF ANALYSIS 


A. THE MODEL 

The aerodynamic model used in this analysis is the same 
as used by DeLaurier and Hui (Ref. 5], axcept that it is ap- 
plied for the lateral case. The assumptions are as follows: 


1) the vehicle is perfectly rigid, flying at a reference 
velocity Us through a constant density op 


ii) the motions are described by the lateral case only 
iii) the control provided by rudder deflection is linearly 
proportional to the yaw angle (wv), that is, 
ACyg 7 kQY 
This last assumption is the case for a helmsman using greater 
control the farther off heading he is perturbed, an assump- 
tion in keeping with operational practice. 
The turbulence model used is the von Karman spectra 
described in Chapter II. It is assumed that the turbulence 
| is composed 2f horizontal gusts only--either Ug or Vg de= 


pending on airship heading relative to the mean wind direc- 


tion. (vg is shown in the analysis.) 


1. Forces from Turbulence Components 


From Jones and DeLaurier [Ref. 6], the normal force 


on a hull segment of length dé is (see figure 6) 


F, = d0U,7(K sin(26) cos| 3] a df + (Cq.),, sing sin|| 2:.dé] 


(8) 


where: K is the hull potential cross-flow factor [Ref. 6] 
6 igs the angle between the hull centerline and Uy 
r is the radius of the hull segment. 

Differentiating with respect to 9 to obtain a perturbation 


equation gives 
ar, = $00," [[- $ sin (ze) sin(3| + 2 x cos(28)cos [5] |aegh ac + 
(Cq.) , (cosdsin| | + sindcos|9|)de ar as] 


By limiting the analysis to the lateral aerodynamic case 
only, (6) becomes (8), the sideslip angle, and Fy becomes 
Yp,r the hull sideforce. Further, if (8,), the undisturbed 
value of sideslip, is assumed to be zero--the usual case-~ 


the above expression becomes 


8 
(4%) 7 x00, KF [2cos(28,)e08(=2|] a0, a 
2. dA 


= oUL K gF dé dé 


Finally, for small values of Var 


Vv 
dg = ma: § 
UN 
and 
2..dA 
(G¥giq = US Kae SE Ct (9) 


29 


The stabilizer forces are given by 


ve (Vos 
(gg * Omgm Sy (Chg) g HUQ) Ng U, 0) 


where: H(k,) is the generalized Sears function as given 
by Filotas [Ref. 16] 


weg ‘ 
kK, = m0, , the "reduced frequency" of the fin 
. (Vols is the gust velocity at the fin mid-chord 


The propellers used to drive the airship produce a 
side force when acted by the turbulence [Ref. 17]. Each 
thruster contribution adds to the total force and moment 
produced, and can be described, for the jth thruster-rotor 


combination by the following: 


U 2 


(Vv ) ' i 
° oy 
i (‘g) 2, = “ep ~Z- Sty ya ty eee tL 1) t 
i 
a | 
(gla, * (Eg)gy (on > Bry) (12) 


Equations (12) and (13) assutne that the rotors are arranged 
symmetrically about the x-2 plane, so that moments due to 


rotor offset in the y-direction cancel out. 


Oe Bie att T Ls nati, . 
2 PLT = 
i tenloani ts at le aes Shek Sl 
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2. Aerodynamic Forces and Moments Due to Airship Motion 


2,qdA,-V av . 
(AY), = 00, ge as ye + oa(k, 4G) + one, Jae (14) 


where Ko, ky are the horizontal and longitudinal apparent 


mass coefficients respectively and 


v(E) sve r{(l, - E)] 


For the fins we have 


(v.) u. s [c z fc yz (c — (15) 
2-9 "0 = + + 
e “2. | ¥B) 5 U5 Yr), 20, p | 
where as ae r(lon - 1,) 
Ve SNe r (lon - 1,) 
(i) ome tS Po he (16) 
Hi he, Oi cat ac Gr 
(Ny) gs * FU Ss, é5{Cne], 20. (17) 


The superscript ac indicates the quantity in parentheses is 
taken about the fin aerodynamic center. 


Thruster forces and moments are given by: 


Uo! (Vn) 4 
Wiylag hee Pay eye a oy, 138) 


Ww 


ydog * hyglgy (Bom - Bry) (19) 


(wp * dy om ~ Bay) (20) 


where 


(Vply =v" rl. = lp,) - p(ho hr) 


3. Inertial Reaction of Airship to Aerodynamic Forces 


and Moments 

The "forces" covered here are those arising as reac- 
tions to airship motion. They are, in general, the negative 
of the forces and moments that give rise to the indicated 
airship translational and angular velocities and accelera- 
tions so as to produce a state of dynamic equilibrium 
(F ~ ma = 0). 


For the hull: 
(4¥ = -¥' (8) (ad, (21) 


where y' =v + r(lon - &) + U5" 


(dL), " "al, + Alor (22) 

(AN), = 7dil,F + aly p (23) 
For the empennage: 

6: es dara. (24) 

(Nady * 7(22) 48 + (TLL) 9P (25) 

(Ls (Tyee) gP = (Toy) gt (26) 
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where aly? ato. and a are the moments and product of 


inertia respectively of the differential element under con- 


sideration, including all structure, air and gas contained 


in the airship. 


4. Bouyancy and Control Terma 


‘ Referring to figure 7, the force due to bouyancy is 
given by: 

a (Yo)_ © ~(gdm - ogAdg) sing cosa, 

sO @? = roll angle 

. a, ™ Steady state angle of attack 


t Differentiating to obtain a perturbation equation gives 


(d¥,) = (ogAdé = gdm) cos}, cosa, ad 


and letting  ™ Q results in 


(d¥ oy ™ [(ogAdé = gdm] cosa, dd (27) 


Control force is assumed to come from rudder deflec- 


tion, and acts through the fin aerodynamic center. 


AC, 


Yo 2 Ka 


2 
Uy 
(Yds 79 "Sy ky 


5. Shear Force, Bending and Twisting Moment 


The hull's shear force at station (1) is obtained by 


(28) 


gumming the sideforce values from the nose, up to (1): 


1 a 
$(l) = [en + ji 8¥a) ay + (Yy) py] (29) 


where (d¥), - (YQ) + (dy + (dY) ny + (d¥y) 


and (a) is the number of rotors forward of station (1). 
Likewise, the bending moment at (1), measured along 


the centerline, is 


1 1 a 
BM(1) = | (Lomrs) (a¥), + | (andy + 3 (bedpy) (%g) gt YW) pd 
i 4 jul J j j 


(30) 
Finally, the twisting moment at station (1) is 
l 1 1 
TM(1) = -| Ha (6) (dY) + | (Alin) » + { (dhig) h 
0 0) 0 
a 
gras) C(Yg) a, + (Yy) ay] (31) 


The ¢carm (aL ia the torque contribution due to the center 


mg) h 
of gravity being offset from the central axis (see figure 7). 


It ig calculated from 


(dog) h = ae dm 6 cosa, 


B. FLIGHT DYNAMICS 


1. Dynamic Stability 
The equations for Lateral Dynamic Stability are taken 


from DeLaurier et. al. [{Ref. 18] and given below: 


CYaero * es - (B - ig) cosa,¢ = 2u(D8 + £) (32) 
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P AC as eu: uchena = 2.8 coma,o = I,,D@ - 1,208 (33) 
7 . AC cece FOCLy.* 2,8 cosa,¢ = Iy,,DB - bf (34) 
Dd =f + 2 tana, (35) 
; DY # seca, (36) 
OCyaero ™ Cy 38 + Cy 408 + Cy, + Cy,Df + Cyn + Cy DP (37) 
4 ACn aero = Cng8 + Cn gd6 + Cnt + Cy pb? + Cn + Cng>P (38) 
| Siaero 7 Clge + C1go8 + Cip# + C1,D2 + C16 + C1508 (39) 
dCy, * Ky (40) 
ee 
cn, * Og = Fon! Cy, ™ Og 7 Fem! KY (41) 
a] G 
| | acy * Poem acy. ™ Coe k GY (42) 
. 
; In this analysis 3 = — = %, thus equations (32) to 
| (39) become ° 
; ACysaro * ACy, ~ (B = fig) cosa, = 2u(De + 2) (43) 
A SCneoro * ACh, 28 cosa, = TDF - 1 DP (44) 
ao 
: | os ar ae 5 Pg 2B cosa, = Ty,DB - Ty,D2 (45) 
po = p + @ tana, (46) 
| Dy = 2 seca (47) 


° 
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ACyaero 
ACn aero 


SC laers 


where: 


[Cc] = 


“ Cy 3? + CyaD0 + Cy .f + Cy,De + Cy, + Cy 506 (48) 


"Cy? + Cygd0 + Caypf + CngDe + Cy, + CnsDB (49) 


= ora + Ci gd + C1 f + C1,be + C156 + C1p>8 (50) 


vis 
ris 
p is 
p is 
y is 


The matrices 


These equations are linear, and along with aquations (40) 


through (42) can be written in matrix form aa: 


> 
oO 
<> 


re> 
Oo 
MH 


- (D] Dp = 0 


€> &> Wh 
Oo 
&> 


translation in the y-direction 


the rotation rate about the z-axis 
the rotation rate about the x-axis 
the roll angle 
the yaw angle 
[C] and [5] are given by 
(Cy .-2u) Cy (mg-B) coc, -k, 
A A (dss lem) 
Cn, Cnp ~X,B cosa, SS aeenG 
a 4 +h 
Clr Slp “2B cosa. EM 
ro) Cc 
tana, l 0 0) 
geCa, 0 ) 0 | 
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[D] = | (2u-Cys) =Cy* “Cy. 0 ) 


B r P 
, “Cha (I27Cny) - (I. 2*Cny) 0 0 
“C14 (Tyo *Cis) (Lyese"C 1g) 0 Q 
7 0 0 0 i 0 
i. 0 0 0 0 1 B 
To solve this system of equations, assume a solution of the 
form 
ee fg ghee. Sw He toe 
| (52) 
43 ee ae y= va ict 
which, when substituted into equation (51) bacomes 
¢ 
R 
i [¢ = ispD] P| = 0 (53) 
| 
| 
‘a ] y 
| where ¢ is the non-dimensional stability root. This is an 
i | eigenvalue problem. A computer golution is performed to find 
a the eigenvalues (stability roots) and eigenvectors (model 
fe. vectors) for the control gains considered. For the subse- 


quent load-response analysis, dynamical: y stable cases must 


be chosen. 
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2. Turbulence Forcing Functions 


As mentioned above, the model chogan for the turbu~ 


lence is sinusoidal with Gaussian statistics, in particular, 


v_(é) casa 
+ = [ exp|iut - ive 7 o (54) 
ro) Q 


By using this expression in equation (9), integrating from 
Se Otogw ly {the hull/fin intersection) and adding the 
contributions of the fins and thruster~rotor combinations, 
the complete turbulence forcing functions for the airship 


can be found. That is, 


a 
L(g) ay + (Yg), (55) 


dy 
Ys | d(Yg)y + 
0 2” 


| 


Likewise, the yawing moment about the nose is 
Nanose ™ +| eee 2 On) Ya) 2] + 1,09), 
Q 


Som * “gnose ~ bem¥g (56) 


and rolling moment is given by 


L. # ~h_¥ (57) 


These may be non-dimensionalized according to 


2 y 2 
ie] ; QO. at 
Yj ey (gig) = 2S (61,6) (58) 
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and the non-dimensional equations can be expressed as 


Gy = Gy G, = Gi,Y as Gny¥ (59) 

where y = [fexp(iwt) = Texp (ike) (60) 
k = ow 
2U 5 


Gy. G1, and Gn, are the turbulent forcing functions for the 
vehicle. 


3. Motion Response Transfer Functions 
Using the functions of equation (59) as forcing 
functions on the right-hand side of equation (53), and 


dividing through by ; gives 


[C = ikD] Fvyr 


G 
y+ 
R/T G 
ny 
p/t = G (61 
ly ) 
o/T ) 


wr | 0 


Solution of this expression for specific geaueed 
frequencies (k), or spectral wave numbers ({), and fixed 
gtable control gains (Ko), allows calculation of the ex- 
pressions necessary for the solution of distributed force 


loadings and moments, by means of the following expressions. 
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‘ Vv a at CG Pp ‘ & 
¢ = a = = exp (iké), ue = exp (ike) 
i = 8 exp(iké), 6 = ¢ exp(ikt) (62) 


= » exp (ike) 
C. LOAD RESPONSE TRANSFER FUNCTIONS 
! Once the motion response of the airship is known, the 
load response transfer functions can be calculated. 
1, Turbulence Loading 
These may be obtained by substituting equation 54 


into equations (9) through (13) and dividing by y (equation if 


| (60)). This gives, for example, for equation (9): 


-_ | (ay Yh ast dA 
Se ht so K SE exp (-iNF cosa.) dé (63) 
r 2 ee 
where 2 = i = 2k 
) c 
Table II gives the complete list of load response transfer ¥ 
functions. ; 
2. Motion Response Loading | 
The aerodynamic-reaction loading may be obtained by 
replacing the motion variables, [+], el, etc., in equa i: 
fs) ° yo 


tions (14) through (20) with the corresponding motion- 
response transfer functions in equations (62). For example, i 


equation (14) becomes 
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i (ay,) U - 6 

a woh o.,dA|V 2 _ +\R 

e pen Bly = ‘tom nF | 

' p 2 c 

i 2 E 2 8 2h 

7 + pU, A{in tS F (tom7 8) | 2 seem k, }}dé 

ye 

(64) 
i: The Inertial=-Reaction, and Bouyancy loading transfer func- 
( tions are similarly obtained from aquations (21) through 

my 

I | (27). For axample, equation (21) becomes 

e > 

t (ay) 2030s : A 

bh m (o} R ’ 2 Vv 2 R 

Eb « (- 2 F + isu, [- - EZ FUtgg 78) | bam (65) 


Coe ot beet 


and equation (27) becomes 


FR 


(daY,,) 
Boh {pgAdg - gdm) cosa, - (66) 


Se A 


; 3. Shear Force, Bending and Twisting Moment Transfer 


MES Pe 


Functions 
The shear force loading transfer function is ob= 


tained using equation (29); 


1 
S(2) dy , 2 (%q)4 
0 jai 
where: 
(¥oln (8% )n | (Xan (En) 


a a aed, wih 
r T r r 


v 
t 
> 
i 


| 
‘| 


BM(1l_) (Y_) (Y,.) 
cm’ gs S ws c's 
aoe + (onde 7 1 2s £ ah 


The bending-moment transfer function comes from aquation (30) 


ZL 1 a (Y_) (Y.4) 
(dy) (an_) qT wT. 
BM(1l) . | (ie ey hy | MR y y (l-1p,) eee thes Cree 
E om r r iat J r r 
0 o - J 
(68) 
and finally the twisting-moment. transfer function, from 
equation (31), is: 
1 1 1 
(dL_) (dL) 
T(t) 7 -| hom (6) (a¥) + | moh , | mg/h 
0 e. & : 
(Y_) (v3) 
a g T wT 
-5 ee + a (69) 
j=1 r r 


An important check on the analytical model is ob- 
tained by ensuring the net Shear Force, and Bending and 


Twisting Moments equal zero. 


S(Uy Was ws | ms | os Ly 


r Tr E r r 


(70) 


The moments are evaluated to the empennage assembly's mass 


center, gt 8° that 


tes) 


c r [ 


——— = 0 (71) 


T (1) (L,,) (L) (Y.) (Y.) (Y.) 
om s , WS, M8 5. (pn ae gs ,0Wws, 08! . 
r r r ig r r r 
(72) 
These equations may be non-dimensionalized as follows: 
1 - 
aoe C. (1) 
k 3 2 $(1 
. = (73) 
: r pus | 
fe) 
Coy (1) oe BM (1) (74) 
iN pu," Se 
Conny (2) | 
r ou, sc 


D. RESPONSE TO ATMOSPHERIC TURBULENCE 

Once the force and moment transfer function coefficients 
are known, the turbulence statistics can be applied to ob- 
tain astimates of airship Lifetime and failure probability. 
When dealing with the lateral aerodynamic case, two different 
Spectra must be considered--$,, and Po9° ‘This analysis will 
use the von Karman spectra as given by aguations (7) and (6) 
respectively. As explained in Chapter II, -the first is used 


when the flight direction is perpendicular to the mean wind, 


and the second when the direction is parallel. 
i 1. Root-Mean-Square Responses 
The root-mean-square response to turbulence of any 


| 
be ; system parameter can be obtained using its transfer function 
y 
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‘ 
sana ees see mes Vs ts ts mae saanonrmeneen armed ini 
TTS BAST LOREM « sib Weaebit tgs ord oe + svat ert natn RE A «cote er 


if multiplied by the spectrum (provided, of course, the spec- 


trum is Gaussian). For example, the rms shear force coef- 
= ficient is | 
a Mefes, eee 1/2 “J 
i (Ca) eng 2 Cy Pia 
cor |) ene an (76) 
2 a 
Cc U. of 


| Response to various conditions can be evaluated using appro- 

! priate values for o and L in the equations for Saye and 
choosing either Pay or 955 according to desired flight 
direction. 


i 2. Mission Analysis Method 


a The mission analysis mathod is a technique for esti- 


sane cl ca 


| mating a flight vehicla's probable lifetime. The method is 7 
based on the probability distribution of encountering tur-~- 


bulence in rapresentative flight operations (Refs. 7 and 15). 


mney cg ae 


{ 
| It is assumed the total flight is a sum of Gaussian patches 
| {Ref. 2]. The formula for calculating the number of ex- 


ceedences is 


-|x-x. | -|x-x_0 | 

a N(x) = JtN_|p,exp ne ee | + piexp ref (77) 
q ojrt b,k : bk 

3 1 2 
x where: x * maximum structural value of bending, moment 
it 


coefficient at a given station, 
x = value of x in one-g level flight, 


N(x) * average number of exceedences of the indicated 


q | value of x per unit time 


Ny = number of zero crossings of x per unit time 

: s . UCase) “Year 

S £ = fraction of time in mission segment 
i 


J Py, Pp * probability values from Table I 


bi, bo = intensity levels from Table I 
also 
; : (Bi) ag] 
m 
yl o tn (BM) 9 | 
where 
pu, 786 2 Coy 
(BM) us = F oe 5 de (78) 
’ 0 
pu, 788} 7/ |Con 
(BM) os = : ia Qs y AL (79) 
0 


{ 
| The probable lifetime is then [N(x)]72 
| 


Shear force and twisting moment are analyzed in the 
game fashion using the appropriate transfer functions. 
3. Other Methods 
DeLaurier and Hui [Ref. 5] also include Failure- 
Probability analysis [Ref. 3] and "Mil-Spec Storm" analysis 


in their paper. There are others in existence, such as the 


Design Envelope Analysis [Ref. 7] that have been used for 
HTA flight, and are well documented. Because of this, they 


| will not be included here, as techniques for using them are 
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the same once the transfer functions of response are known. 


The exact method used ig up to the designer, based on his 


needs. 
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IV. NUMERICAL EXAMPLE 


In order to illustrate the lateral aarodynamic case 
developed in Chapter III, an example using the USS AKRON 
(ZR=4) is presented. The flight conditions chosen are: 

U, = 123 ft/sec 


Alt = 1000 ft 
The velocity represents the maximum for the vehicle, and the 
altitude is typical of its operational range. In addition, 
a condition of neutral buoyancy (B-mg = 0.0) was selected. 
The geometry was taken from Freeman [Ref. 19] and the weight 
distribution from Woodward [Ref. 20]. With this information 
available, the inertial properties of the AKRON could be 
calculated using the method of Scholaert and DeLaurier [Ref. 
21] (see Appendix). The values obtained are shown in Table 
Til. The Hull cross-flow and stabilizer efficiency factors 
are calculated using the method given in reference 6, and 
are: 

K = 0.93225 
Ng ™ 9.2600 


The apparent-mass coefficients are from Munk [Ref. 22]. 
The stability derivatives are taken from DeLaurier and 
Schenck [Ref. 18], and shown in Table IV. With these, the 


control gain, and the inertial and geometrical properties, 


aquation 53 can be solved to find the stable roots. This 
was done in reference 18 and the rasults are shown in figure 
8. Mode 4 (indicated in the figure) is characterized by 
roll, yaw, and gideslip of equal magnitudes not unlike the 
dutch roll mode of a fixed wing aircraft. Mode 5 is one of 
equal and opposite @ and ) motions with small » perturka-~ 
tions. Mode 6 is a relatively high-frequency rolling motion, 
little affected by control gain. Modes 1, 2, and 3 refer to 
longitudinal aerodynamic modes [Ref. 18]. From this analy- 
gis, a control gain of 0.2 was found to provide the minimum 
stable condition. 

The forcing functions were next obtained using equations 
§4 through 60, and are plotted in figure 9. The peaks in 
all the curves occur at a wave number of about .008. This 
corresponds to a wavelength equal to the length of the air- 
ship. $1, | is significantly smaller than the others, as 
expectad, due to tha smaller moment arm through which the 
sideforce works in producing roll. 

With the turbulence forcing functions, equation 61 was 
solved to obtain tha motion response transfer functions. 
Thase are shown in figures 10 through 14. Control gains of 
0.2, 1.0, and 2.0 were used to illustrate the effect of its 
variation. The most significant feature of these responses 
is the peak at a wave number of .008. This corresponds, as 
in the forcing functions, to a condition where the spectral 


component wavelength exactly equals the airship length. This 
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is the result predicted by Calligeros and McDavitt [Ref. 4] 
for the longitudinal case. DeLaurier and Hui [Ref. 5] also 
obtained this result, but only for cases of higher control 
gain. For the lateral case, the control gain does not sig- 
nificantly change this peak, although, for yaw and yaw rate 
(y and y), and to a lesser extent roll (¢), the response at 
lower wave numbers is reduced. 

Finally, the load response transfer functions were cal- 
culated. The results are shown in figures 15 through 17 for 
a wave number of .009, and a control gain of 0.2. Complete 
results are given in the appendix. For the most part, the 
results yield no surprises. The magnitudes of the load 
response follow the general trend of the combined motion 
responses, thug the peak loads occur at wave numbers near 
.008. The location along the axis of the peak load varies 
as the magnitude of the motions increases, shifting aft in 
the case of shear and twisting moment, and to the center for 
bending moment. Again, the lack of significant change with 


control gain is apparent. 


ee 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The analysis presented is an extension of the work by 
DeLaurier and Hui [Ref. 5], and is subject to the same re-~ 
strictions. That is, it is limited to small perturbations 
in order to allow a linear analysis usable with power spec- 
tral methods, and its ability to make precise predictions of 
the loading when used with one of the methods that accounts 
for severe turbulence is questionable. Nonetheless, it is 
a valuable tool in understanding the response to an initial 
disturbance, and when employed as the aerodynamic input to 
the various statistical methods discussed earlier, it can 
yield important design and operational insight. 

The limited effectiveness of the simple control model 
employed, which ig typical of someone cuing his response to 
a compass, was demonstrated. It is suggested that an exami- 
nation of the effects of roll control and yaw rate feedback 
be made. This would allow a decision as tc the feasibility 
of using control to provide gust alleviation. As discussed 
by DeLaurier and Hui, control gain made a large difference 
in the expected lifetime of an airship when considering only 
longitudinal aerodynamics. Undoubtedly, for the lateral 
case, aven the yaw control used in this analysis will con- 


tribute to increased survivability due to the reduction of 


loads at low wave numbers. Thus, the next step is to employ 


the statistical methods to discover how mich change is 
realized. 

The cuse of combined longitudinal and lateral motion 
needs to be studied. No aircraft ever built Fas ever 
managed to fly through turbulence tha. is strictly one 
dimensional, as is asswned for this analvsis. Coupling the 
two cases would give a much better idaa of the true action 
of an airship in turbulence. 

Finally, some means must be found to astablish the vera~ 
city of this model, as well as that for the longitudinal 
case. To the author's knowledga, no investigation of the 
actual response of an airship to conditions of known turbu- 
lence has ever been made. this is, of course, a difficult 
project, considering the limited number of airships cur- 
rently available if full scale tests are to be carried out. 
Wind tunnel investigations, made in the various oscillating 
flow tunnels available, would be helpful. Until some tests 
are done, howaver, this type of analysis must be considered 
only for its qualitative aspects as opposed to its quantita- 


tive predictions. 
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COMPUTER OUTPUT FOR THE NUMERICAL EXAMPLE 
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TABLE IIt 
GEOMETRICAL AND INERTIAL PROPERTIES 


OF THE USS AKRON (ZR=4) 


total volume = 7,382,400 £t° 
¢ = 785.0 ft 
s = 37,914 £t* 


1am m 364.24 fit 


Dom = -37.66 ft 


L, = 363.01 ft 


mass = 17,039 slugs 
buoyancy =# 548,642 lb 


Tee % 38,685,100 slug-£t? 


Tez * 471,799,000 siug-£t? 


Ip, * 102,129,000 slug-ft? 
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* TABLE IV 


STABILITY DERIVATIVES OF THE USS AKRON (2R=4) 


neutral buoyancy, U. = 123 ft/sec, ALT = 1000 ft 


° 
i Cy, = -0.7224 Cy, = -0.9863 
Cy, = 70.3418 Cy, = 0.0586 
| Cy, = -0.0648 Cy, = -0.0913 
Cng = 70.1710 Chg = 0.0293 
Ch. = 0.2352 Che ™ -0,0991 
Cnpy = 0.0150 Chg * 0.0028 
C1, = 70.0322 Cy = 70.0456 
C1. = -0.0153 Cy, = 0.0028 


| Cy. = -0.0066 Cy, = -0.0042 
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' Figure 1. The (1~Cosine} Gust Shape 
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Figure 4. Typical Power Spectra of Vertical Gust Velocity 
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, Figure 5. Measured and Fitted von Karman Spectra 
a ; of Vertical Gust Velocity from Severe Storm 
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Figure 7. Schematic of Buoyancy Forces and Moments 
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Figure 8. Lateral Root-Locus of the USS AKRON (ZR-4) 


138 


suoT}OUN, HuToi0g asoustnqing, °6 eAznbta 


(33/suetpe1r) rA5qumN SAP 


ot 7-01 5-01 ol OT 


rd 
‘ 


Oot 


guotqzound butozo0d souetnqang, 


139 


esuodssy dtis eprs ‘oT eanbra 


(33/suetpez) zaequnyg oaem 
-ot 7-01 ,-ol pot gol 
———=., 19° 


0°z = (4 —--— w 
O°l = 24----- > 
zZ°0 = °% l' wo 
| ol 
p- 
ue) 
o = 
3 wt 
it] 
_ ie 
. T 
A<> 
oT 


oo ghee ee cae ge oT ae 


[4 


esuodsoy MAPA “TT eanbta 


(ay /suetper) aequmN ocAeM 
_Ot OT 


red 
° 


ot 


‘asuodsey MB 


E 


141 


esuodsay 3o3ey AMPK “ZT eanbtg 


{33 /SueTpert) ASquNN eAeM 
oT oT OT y-01 ¢-0T 


ed 
® 


f‘asuodsey ater MeZ 


OT 


142 


esuodsey IfTow “¢1T eanbra 


(43 /SueTper) zASqUMN SAeM 


7-01 z_ol ¢—OT ,-ol g-OT 
TO 
wa 
o 
T° at 
we 
wD 
oo 
8 
tc 
g 
Can] 
a 
® = 


ote om gnaeee 


asucdsoy o2eu TTOWw “FT oeanbta 


(az/suetpez) tequmN eAeM 
ot z_ot ¢-ot 4~ot Ol 


el 


144 


f-~] 
rm 
‘gsuodsey ave TTOu 


i 


oot 


00s 


OOL 


queTSOTZJ2ZO0D soAO, aAveYS 


(43) 3 
009 00s oor 


“GT oan6Bta 


coe 00¢c 


600° 
z°0 
Des/FF EZT 
“33 000T = 


Oot 


i) 
re 


wl 
=~ 


° 
™N 


Se 


‘4QUOTOTS FOOD VOAOY rzveYS 


"3 


145 


qUueToTJye0D yueMCH Gurpueq “9{[ eanSTi 


(az) 3 
008 OOL 609 90S 00F 00t 00z 00T 1) 
7) 
w 
oO 
3 
is 
ft “ 
n 
K4 
oO 
| 
3 
4 t 
Q ‘ 
fe) = 
th et 
th 
£ P- 
Q 
600° = 3 4 
] ro] 
e0-= A c 
oO 
oesf/AF €ZT = A ce 
Q 
“233 OOOT = *3TV — 
G 


Pe RRB ROR Ra a ee ee ee ee one ee. 


quaToTyjooD USTIOM HUTASTML “LT eanbta 


(93) 3 
008 002 009 06S 00V O0L god oot 


44quaeTotTyyJe0D yZusWOP HuUTIaTAL 


S°0 


147 


600° = B 
zo= 4 
Ses/f3 €ZI = A 
“33 900T = “FIV 0°T 


LIST OF REFERENCES 


ve Burgess, C. P., Airship Design, pp. 99-100, Ronald 
l Press, 1927. 


2. Houbolt, J, C., Steiner, R., and Pratt, K. G., Dynamic 
Response of Airplanes to Atmospheric Turbulence In- 
sliding Flight ve on input and Rasponse, NASA TR-R- 
199, ry 

3% Press, H., and Meadows, M. T., A Reevaluation of Gust- 
Load Statistics for Applications in Spectral Calcula- 
tilons, NACA TN 3540, EEE 

4. Calligeros, J. M., and McDavitt, P. W., Response and 
Loads on Airships due to Discrete and Random Gusts, MIT 
Aeroelastic and structures Research hab. Tech. Rept. 
72-1, Cambridge, Massachusetts, February 1958. 

| ' =a DeLaurier, J. D., and Hui, K. Cc. K., Airship Survivabilit 

in Atmospheric Turbulence, American Institute of Aero- 

nautica and Astronautics Paper No. 81-1323, July 1981, 

; 6. Jones, S. P., and DeLaurier, J. D., Aerodynamic Estima- 

tion Techniques for Aerostats and Airships, ALAA Paper 

No. 81-1339, July 1981. 

| 7. | Stkin, B., The Turbulent Wind and Its Effect on Flight, 

| UTIAS Review No. 44, University of Toronto institute 

| for Aerospace Studies, August 1980. 

8. Batchelor, G. K., Theory of Homogeneous Turbulence, 
Cambridge University Press, Cambridge, Massachusetts, 


1953. 

9. Dobrolenskiy, Yu. P., Flight Dynamics in Moving Air, 
NASA TY F-600, July 1971. 

10. Etkin, B., Dynamics of Atmospheric Flight, Wiley, 1972. 


a . 11. Ferzinger, J. H., Large-Eddy Simulation, Lecture pre~ 
| sented at the AIAA Professional Study Seminar on Tur- 
bulence Modeling, Palo Alto, California, 20-21 June 1981. 


12. Ribner, H. S., Spectral Theery of Buffeting and Gust 
Response; Unification and Extension, Journal of Aero. 


| Scl., Vol. 23, No. 12, 1956. 


148 


13. 


14. 


1 Gs 


16. 


17. 


18. 


19. 


20. 


21, 


22. 


Air Force Flight Dynamics Laboratory, An Exposition on 
Aircraft Response to Atmospheric Turbulence Using Power 
Spectral Density Analysis Techniques, Technical Report 
AFFDL-TR-7 0-16 7 ‘May dd oes 


Houbolt, J. C., “Atnospheric Turbulence," AIAA Journal, 
Vol. 11, No. 4, April 1973. 


National Aeronautics and Space Administration, Terres- 
trial Environment (Climatic) Criteria Guidelines for 
Use in Aerospace Vehicle Deve opments Revision, 
edited by Kaufman, J. W., NASA TM 78118, November 1977. 


Filotas, L. T., “Approximate Transfer Functions for 
Large Aspect Ratio Wings in Turbulent Flow," AIAA J. of 
Aircraft, Vol. 8, No. 6, June 1971. 


Ribner, H. S., Propellers in Yaw, NACA Report No. 820, 
1945. 


DeLaurier, J. D. and Schenck, DP. M., Airship Dynamic 
Stability, ALAA Paper No. 79-1591, July ror 
Freeman, Hugh B., Force Measurements on a 1/40=Scale 


Model of the U.S. Alrship ‘Akron’, NACA Report No. 432, 
1932. 


Woodward, Donald E., Private Communication, Association 
of Balloon and Airship Constructors, 1981. 


Scholaert, H. and DeLaurier, J. D., Private Communica- 


tion, University of Toronto Institute for Aerospace 
Studies, 1981. 


Munk, Max M., The perogppamis Forces on Airship Hulls, 
MACA Report No. ’ . 


‘s INITIAL DISTRIBUTION LIST 


No. Copies 
1. Defense Technical Information Center 2 
Cameron Station 
Alexandria, Virginia 22314 
2. Library, Code 0142 2 
Naval Postgraduate School 
Monterey, California 393940 
3. Department Chairman, Code 67 an 
Department of Aeronautics 
Naval Postgraduate School 
Monterey, California 93940 
4, Professor Donald M. Layton, Code 67Ln 10 
Department of Aeronautics 
Naval Postyraduate School 
Monterey, California 93940 
5. LT John J. Wrobleski, Jr. 5 


USS DWIGHT D. BKISENHOWER (CVN-69) 
FPO New York 09532 


